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This  study  uses  the  China  Coastal  Waters  and  Adjacent  Seas  Reanalysis  (CORA)  data  to  investigate  the  interannual  and  decadal 
variability  of  the  Yellow  Sea  cold  water  mass  (YSCWM)  and  its  relationship  to  climate  indices  including  the  Arctic  Oscillation  (AO), 
El  Nino/Southern  Oscillation  (ENSO),  and  Pacific  Decadal  Oscillation  (PDO).  On  the  interannual  timescale,  the  strong  correlation 
between  winter  PDO  and  the  YSCWM  may  indicate  the  dominant  effect  of  winter  PDO  on  the  YSCWM  through  the  modulation 
of  local  heat  flux  and  wind  stress.  It  is  also  found  that  the  local  wind  stress  and  heat  flux  in  summer  have  little  impact  on  the 
interannual  variability  of  the  YSCWM.  On  the  decadal  time  scale,  the  YSCWM  is  associated  with  winter  AO  and  winter  PDO. 
Winter  AO  mainly  controls  local  heat  flux,  modulating  the  decadal  variability  of  the  YSCWM.  In  contrast,  winter  PDO  is  strongly 
connected  with  winter  heat  flux  and  wind  stress  to  modulate  the  decadal  variability  of  the  YSCWM.  In  summer,  for  three  climate 
factors,  ENSO  is  the  dominant  factor  controlling  the  decadal  variability  of  the  YSCWM. 


1.  Introduction 

The  Yellow  Sea  is  a  semienclosed  marginal  sea  of  the  Pacific 
Ocean,  bordered  by  the  Chinese  mainland  to  the  west,  the 
Korean  Peninsula  to  the  east,  and  the  Bohai  Bay  to  the 
north.  Several  rivers  drain  into  the  Yellow  Sea  and  provide 
a  lot  of  nutrients.  The  Yellow  Sea  covers  an  area  of  roughly 
400,000  km^  and  it  is  mostly  shallower  than  80  m  (Figure  1). 
The  central  trough  extends  from  southeast  to  north  with  an 
average  depth  of  44  m  [1]. 

The  seasonal  variation  is  evident  in  the  circulation  and 
hydrographical  characteristics  of  the  Yellow  Sea.  In  winter, 
strong  vertical  mixing  due  to  strong  winds  and  surface 
heat  fluxes  produce  a  uniform  and  cold  vertical  temperature 
profile.  In  summer,  the  solar  radiation  increases  and  rapidly 
warms  the  upper  layer  of  the  Yellow  Sea.  Therefore,  the 
seasonal  thermocline  rapidly  forms  and  locates  at  10-30  m 


depth.  Due  to  the  combined  topographic,  dynamic,  and 
thermodynamic  effects,  there  is  a  remnant  of  cold  water 
(<10°C)  under  the  seasonal  thermocline  in  the  central  trough 
of  the  Yellow  Sea  which  is  referred  as  Yellow  Sea  Cold 
Water  Mass  (YSCWM)  and  is  present  throughout  the  whole 
summer  season  occupying  almost  30%  of  the  Yellow  Sea  area, 
thus  becoming  one  of  the  most  important  characteristics  of 
the  Yellow  Sea.  The  YSCWM  has  an  important  effect  on  the 
hydrographic  features  and  the  phytoplankton  biomass  and 
production  in  the  Yellow  Sea  [2-5]. 

He  et  al.  [6]  first  studied  the  formation  and  properties 
of  the  YSCWM  and  identified  that  the  water  mass  locally 
formed  during  the  previous  winter  by  the  sea  surface  cooling 
and  strong  vertical  mixing.  Su  [7]  found  that  the  bottom 
cold  water  in  northern  Yellow  Sea  is  influenced  by  coastal 
water  and  then  moves  southward  in  winter,  mixing  with  high- 
salinity  water  from  the  open  ocean.  The  bottom  cold  water  in 
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Figure  1:  Bottom  topography  of  the  Yellow  Sea  (m). 


northern  Yellow  Sea  remains  at  the  bottom  in  summer,  and  it 
is  influenced  by  the  Yellow  Sea  Coastal  Water  [8].  Since  1990, 
numerical  models  are  used  to  investigate  the  mechanism 
of  the  YSCWM  formation  and  its  associated  circulation 
[9-13].  Su  and  Huang  [9],  using  qualitative  analysis  and 
numerical  simulations,  obtained  that  the  vertical  circulation 
is  characterized  by  a  double  cell  structure.  Miu  et  al.  [10, 
11]  used  a  model  to  explain  the  formation  of  the  YSCWM 
and  basic  characteristics  of  density  circulation  based  on  the 
fluid  dynamics  and  the  thermodynamics.  Recently,  several 
studies  were  conducted  to  investigate  the  characteristics  of 
the  YSCWM  using  temperature  and  salinity  observation. 
Jiang  et  al.  [14]  studied  its  thermohaline  characteristics  based 
on  the  temperature  and  salinity  data  observed  by  the  State 
Oceanic  Administration  Peoples  Republic  of  China  from 
1976  to  1999.  Hu  and  Wang  [15]  used  empirical  orthogonal 
function  to  investigate  the  interannual  variability  with  tem¬ 
perature  data  collected  in  August  from  1975  to  2003.  Although 
previous  studies  investigated  the  interannual  variability  of 
the  YSCWM  through  some  frequently  measured  temperature 
sections,  due  to  the  lack  of  a  long  and  regional  observational 
record  they  could  not  describe  the  interannual  and  decadal 
variability  at  a  regional  scale.  However,  ocean  reanalysis  data 
has  a  better  space-time  coverage,  which  allows  us  to  inves¬ 
tigate  the  whole  features  of  the  YSCWM  and  its  response  to 
climate  variability.  Because  the  YSCWM  is  the  most  conser¬ 
vative  among  water  masses  in  the  Yellow  Sea,  it  may  contain 
clearer  climate  signals  than  any  other  water  masses  in  the  Yel¬ 
low  Sea.  These  long-term  signals  are  important  to  understand 
the  climatological  evolution  of  the  Yellow  Sea  [16]. 

Therefore,  this  study  aims  to  characterize  the  low- 
frequency  variability  of  the  YSCWM  and  it  possible  rela¬ 
tionship  with  three  climate  indices:  the  Arctic  Oscillation 
(AO),  El  Nino/Southern  Oscillation  (ENSO),  and  the  Pacific 
Decadal  Oscillation  (PDO)  using  the  China  Coastal  Waters 


and  Adjacent  Seas  Reanalysis  (CORA)  data.  CORA  assimi¬ 
lated  the  in  situ  observed  temperature  and  salinity  profiles 
collected  by  the  National  Marine  Data  and  Information 
Service  Peoples  Republic  of  China,  World  Ocean  Data  2005 
(WOD05),  Global  Temperature  and  Salinity  Profile  Project 
(GTSPP),  and  Array  for  Real-Time  Geostrophic  Oceanogra¬ 
phy  (Argo)  floats.  It  is  worthy  to  emphasize  that  the  dataset 
from  the  National  Marine  Data  and  Information  Service 
Peoples  Republic  of  China  integrates  almost  ah  the  available 
historical  observed  temperature-salinity  profiles  in  the  Yellow 
Sea,  which  contributes  to  the  high  quality  of  CORA  in  the 
Yellow  Sea.  CORA  has  been  validated  using  independent 
observations  in  previous  studies.  The  results  showed  that 
CORA  can  better  simulate  temperature-salinity  structure, 
mesoscale,  and  large-scale  circulation  in  the  western  north 
Pacific  including  the  Yellow  Sea  [17] .  Eor  further  details  about 
CORA  and  the  assimilated  observations  the  reader  is  referred 
to  [17]. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section  2, 
we  describe  the  data  used  in  this  study.  In  Section  3,  we 
explore  the  seasonal  evolution  and  formation  of  the  YSCWM. 
In  Section  4,  we  investigate  the  interannual  and  decadal  vari¬ 
ability  of  the  YSCWM  and  discuss  the  relationship  between 
the  YSCWM  and  the  three  climate  indices  (AO,  ENSO, 
and  PDO)  considered  in  this  study.  Section  5  gives  a  brief 
summary  and  discussion. 

2.  Data  and  Methods 

The  CORA  covers  a  period  of  51  years  from  January  1958  to 
December  2008.  The  domain  ranges  from  99°E  to  I50°E  and 
from  I0°S  to  52°N,  including  the  Bohai  Sea,  the  Yellow  Sea, 
the  East  China  Sea,  the  South  China  Sea,  and  adjacent  seas. 
The  datasets  consist  of  monthly  mean  fields  with  horizontal 
resolution  of  0.25°  and  35  vertical  depth  levels. 
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Figure  2:  Monthly  mean  climatology  of  bottom  water  temperature  (°C)  in  the  Yellow  Sea.  The  triangles  indicate  the  location  where  the  cross 
sections  are  performed. 


To  examine  what  climate  factors  are  associated  with 
the  YSCWM  variability,  three  climate  indices  are  selected: 
AO  index  [18],  Multivariable  ENSO  index  (MEI)  [19,  20], 
and  EDO  index  [21,  22].  AO  is  the  dominant  pattern  of 
nonseasonal  sea  level  pressure  variability  north  of  20° N.  It  is 
characterized  by  sea  level  pressure  anomalies  with  opposite 
signs  in  the  Arctic  and  37-45°N  [18].  MEI  is  an  average  of 
the  main  ENSO  features  contained  in  six  observed  variables 
over  the  tropical  Pacific.  A  positive  (negative)  value  of  the 
MEI  indicates  the  warm  (cold)  ENSO  phase  [19,  20].  The 
PDO  index  is  defined  as  the  leading  principal  component  of 
monthly  sea  surface  temperature  variability  north  of  20° N  in 
the  North  Pacific,  and  it  is  a  long-lived  El  Nino -like  pattern 
of  the  Pacific  climate  variability  [21,  22]. 

The  monthly  mean  wind  stress  and  net  heat  flux  are 
computed  averaging  the  daily  values  obtained  from  the 
NCEP-NCAR  reanalysis  with  the  resolution  of  2.5°  [23].  The 
period  of  51  years  ranging  from  January  1958  to  December 
2008  is  used  to  investigate  the  mechanisms  of  the  YSCWM 
variability  associated  with  AO,  ENSO,  and  PDO.  To  examine 
the  interannual  and  decadal  variability  of  the  YSCWM, 
all-time  series  used  in  this  study  are  separated  into  two 
components  including  inter  annual  and  decadal  variability. 
The  decadal  component  is  obtained  by  applying  a  7-year 
Gaussian  filter  to  the  original  time  series.  The  interannual 


component  is  obtained  by  subtracting  the  decadal  component 
from  the  original  time  series.  To  extract  the  dominant  mode 
of  heat  flux  and  wind  stress  in  the  Yellow  Seas  region,  an 
empirical  orthogonal  (EOE)  analysis  is  performed  upon  the 
interannual  and  decadal  time  series. 

3.  Seasonal  Evolution  and 
Formation  of  the  YSCWM 

The  YSCWM  can  be  clearly  distinguished  in  the  bottom  of 
the  Yellow  Sea.  The  Yellow  Sea  bottom  water  is  the  most  stable 
level  of  the  YSCWM  [6].  YSCWM  is  traditionally  defined  as 
the  cold  water  lower  than  10°  C  in  the  bottom  of  the  Yellow 
Sea  in  summer  [6].  Thus,  climatological  monthly  averaged 
temperature  of  the  Yellow  Sea  bottom  water  during  1958- 
2008  is  used  to  examine  the  formation  and  evolution  of  the 
YSCWM  (Eigure  2). 

In  winter  (December-Eebruary),  the  temperature  is  grad¬ 
ually  reduced  in  the  Yellow  Sea.  The  cold  water  area  (with 
temperatures  below  10° C)  is  gradually  increased.  In  January, 
there  is  a  tongue  of  warmer  water  (with  temperatures  warmer 
than  10°C)  whose  tip  moves  from  (35°N,  122.5°E)  to  (34°N, 
123°  E)  from  January  to  Eebruary  and  further  moves  south¬ 
eastward  in  spring  (Eigure  2).  The  temperature  in  the  central 
trough  is  higher  than  that  of  the  surrounding  water  in  winter. 
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Figure  3:  Continued. 
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Figure  3:  Vertical  cross  sections  of  monthly  mean  temperature  (°C)  along  36.17°N  and  38.667°N  for  the  locations  indicated  in  Figure  2. 


In  spring  (March-May),  the  Yellow  Sea  water  is  warming 
further,  and  its  frontal  structure  is  changing.  In  March,  the 
temperature  structure  is  similar  to  that  in  February.  The  warm 
water  tough  still  exists  and  then  gradually  disappears  in  April. 
However,  in  May,  there  are  three  cold  water  patches  with 
temperatures  below  7°C  form  around  36°  25'N,  122.5°E,  35° 
50^N,  123.5°E,  and  36°  25^N,  125°E  (Eigure  2).  This  suggests 
that  the  YSCWM  has  started  to  form  in  the  Yellow  Sea. 

In  summer  (June- August),  the  cold  center  gradually 
enhances.  In  June,  the  cold  center  located  at  36°  25^N,  122.5°E 
starts  to  move  eastward,  while  the  other  two  cold  water 
centers  merge  into  one  during  June,  suggesting  that  the  main 
body  of  the  YSCWM  begins  to  form  over  southern  Yellow 
Sea.  The  cold  water  mass  in  the  bottom  of  the  central  trough 
has  a  shape  that  resembles  a  saddle.  In  July,  the  cold  water  that 
has  a  temperature  lower  than  10°  C  covers  a  third  of  the  Yellow 
Seas  total  area.  Compared  with  that  in  the  late  spring,  the 
isotherm  of  10° C  expands  southward  from  36°N  to  35.25°N 
(Eigure  2).  This  indicates  that  the  YSCWM  has  reached  it 
maximum  extension. 

In  autumn  (September-November),  the  temperature  of 
the  coastal  water  drops,  especially  in  northern  Yellow  Sea. 
Although  the  YSCWM  is  still  visible  until  November,  the 
lowest  temperature  at  its  center  is  around  12°  C,  and  its  area 
has  also  reduced.  In  December,  the  YSCWM  has  completely 
disappeared. 


We  have  analyzed  the  formation  and  evolution  of  the 
YSCWM  through  the  horizontal  distribution  of  temperature 
at  the  Yellow  Sea  bottom.  Next,  we  choose  two  vertical 
temperature  sections  of  the  YSCWM  across  southern  and 
northern  Yellow  Sea,  respectively.  The  location  of  the  two 
sections  is  indicated  with  triangles  in  Eigure  2.  We  will  further 
explore  the  formation  and  the  evolution  mechanism  of  the 
YSCWM  by  analyzing  the  vertical  structure  of  temperature 
in  different  seasons. 

Over  northern  Yellow  Sea,  the  cross  section  is  performed 
at  38.66°N  through  the  center  of  the  YSCWM.  As  shown  in 
Eigure  3,  in  winter,  the  water  column  is  completely  uniform. 
In  April,  the  enhanced  solar  radiation,  the  weaker  convection, 
and  mixing  destroy  the  vertical  uniform  distribution  of 
temperature,  and  a  thermocline  is  formed  due  to  water 
stratification.  But  this  thermocline  is  weak  with  a  maximum 
gradient  of  ~0.2°C/m  at  approximately  10  m  depth.  In  May, 
the  thermocline  develops  rapidly,  and  its  intensity  is  greater 
than  0.3° C/m  at  approximately  15  m  depth.  At  the  same  time, 
an  area  with  temperatures  under  8°C  forms  underneath  20  m 
with  minimum  temperatures  under  7°C  near  38.667°N,  123°E. 
This  cold  center  is  also  visible  in  Eigure  2.  The  YSCWM  has 
started  to  form  at  the  bottom  in  May. 

In  early  summer,  the  stratification  at  depths  between  10 
and  20  m  became  evident,  and  the  thermocline  becomes  sta¬ 
ble.  The  temperature  profile  presents  a  three-layer  structure: 
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Figure  4:  Interannual  time  series  of  (a)  winter  AO,  MEI,  and  PDO  indices,  and  A-YSCWM  and  T-YSCWM  in  August  (b)  summer  AO,  ENSO, 
PDO,  and  the  A-YSCWM  and  T-YSCWM  in  August.  Time  series  of  the  A-YSCWM  and  T-YSCWM  are  standardized.  The  A-YSCWM  and 
T-YSCWM  represent  the  area  and  averaged  temperature  of  the  YSCWM,  respectively.  The  “+”  represents  the  El  Nino  events,  and  the 
represents  the  La  Nina  events. 


the  upper  mixed  layer  contains  warm  water;  the  thermocline 
has  sharp  temperature  gradient  (0.2-0.5°C/m);  and  the  lower 
layer  contains  cold  water  which  is  nearly  uniform.  In  mid¬ 
summer,  the  sea  surface  temperature  rises  rapidly,  and  the 
strength  of  thermocline  reaches  its  maximum,  which  is  as 
high  as  1.0° C/m.  At  the  same  time,  the  bottom  water  still  has 
temperatures  below  10° C.  In  autumn,  the  temperature  of  the 
cold  center  rises  from  8°C  to  11°  C  (Figure  2).  This  indicates 
that  the  YSCWM  is  gradually  decaying.  In  November,  the 
thermocline  disappears  completely,  and  the  YSCWM  is  also 
gone. 

We  choose  the  section  at  36.167°N  to  examine  evolution  of 
the  YSCWM  in  southern  Yellow  Sea  (Figure  3).  The  evolution 
process  of  the  thermocline  in  the  southern  Yellow  Sea  is 
similar  to  the  northern  Yellow  Sea,  but  due  to  the  different 
geographical  position,  topography,  and  tidal  mixing  intensity, 
their  vertical  structures  are  different.  In  spring,  as  the  tem¬ 
perature  increases  near  the  coast,  the  thermocline  weakens 
on  both  coasts  of  the  southern  Yellow  Sea,  and  the  area  of  the 
cold  water  increases.  In  May,  with  the  increasing  temperature 
and  enhancing  thermocline,  three-layer  structure  occurs.  The 
westward  moving  of  the  east  cold  center  makes  cold  water 
only  exists  in  the  central  trough.  In  June-August,  the  tem¬ 
perature  of  the  bottom  water  in  the  central  trough  is  less  than 
10° C.  In  October,  the  thermocline  only  exists  at  the  central 
deeper  water  area  than  30  m.  The  intensity  of  thermocline 
also  reduces  significantly.  By  December,  the  thermocline 
completely  disappears,  and  the  YSCWM  is  also  gone. 

There  are  two  cold  centers  below  the  thermocline  located 
at  122. 5°E  and  124.5°E  in  the  southern  Yellow  Sea  in  April, 
which  are  different  at  the  northern  Yellow  Sea  (Figure  3). 
This  difference  is  due  to  the  influence  of  the  Yellow  Sea 
warm  current.  In  winter,  a  warm  tongue  from  west  of  the 
Jeju  Island  stretches  northward  along  the  central  trough 
(Figure  2),  which  causes  the  temperature  over  the  central  part 
to  be  higher  than  that  over  the  two  coasts  of  the  Yellow  Sea. 
For  example,  in  the  cross  section  performed  at  36.167°N,  the 
temperature  in  the  central  part  of  the  southern  Yellow  Sea  is 
7°C  in  February.  However,  the  temperature  on  both  coasts  is 
4~6°C.  In  April,  the  temperature  over  the  central  region  is 
greater  than  8°C,  but  the  temperature  on  both  coasts  is  6~7°C 


(Figure  3).  The  bottom  temperature  in  the  southern  Yellow 
Sea  has  two  cold  centers  on  both  flanks  of  central  trough  from 
January  to  April  (Figure  3). 

4.  Interannual  and  Decadal  Variability 
of  the  YSCWM 

The  YSCWM  is  a  large  water  mass  covering  a  third  of  the 
bottom  layer  at  its  largest  extension.  It  has  relatively  stable 
properties  with  low  temperature  and  high  salinity.  However, 
due  to  the  influence  of  external  forcings  such  as  wind  stress 
and  heat  fluxes,  the  thermohaline  structure  within  the  water 
mass  may  change.  The  YSCWM  is  traditionally  defined  as  the 
cold  water  lower  than  10° C  in  the  bottom  of  the  Yellow  Sea 
in  summer  [6].  To  investigate  the  interannual  and  decadal 
variability  of  the  YSCWM,  we  define  the  YSCWM  as  the 
bottom  water  colder  than  10°  C.  To  examine  the  effects  of 
the  different  seasonal  climate  forcings  on  the  interannual 
and  decadal  variability  of  the  YSCWM,  winter  (summer) 
climate  indices  averaged  from  December  to  February  (July  to 
September)  are  used. 

4.1.  Inter  annual  Variability 

4. 1.1.  Relationship  to  the  Climate  Forcings  in  Winter.  On  the 
interannual  time  scale,  both  winter  AO  and  ENSO  indices  are 
negatively  (positively)  correlated  with  the  area  and  averaged 
temperature  of  the  YSCWM  (hereafter  referred  to  as  A- 
YSCWM  and  T-YSCWM,  resp.).  Although  the  correlation 
coefficients  are  low,  it  is  clearly  shown  that  El  Nino  events 
(e.g.,  1960,  1983,  1992,  and  1997)  correspond  to  smaller  A- 
YSCWM  and  higher  T-YSCWM,  while  La  Nina  (1963, 1974, 
1985)  events  correspond  to  larger  A-YSCWM  and  lower  T- 
YSCWM  (Figure  4(a)  and  Table  1).  The  correlation  between 
the  winter  PDO  index  and  A-YSCWM  (T-YSCWM)  is  pos¬ 
itive  (negative)  and  statistically  significant  at  the  confidence 
level  of  95%.  The  correlation  between  the  winter  PDO  index 
and  the  YSCWM  is  the  strongest  among  all  of  the  three  cli¬ 
matic  indices,  which  may  indicate  that  the  effect  of  winter 
PDO  on  the  YSCWM  is  the  dominant  climate  forcings  (Fig¬ 
ure  4(a)  and  Table  1). 
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Figure  5:  EOFl,  EOF2,  and  EOF3  of  interannual  component  for  winter  heat  flux  ((a),  (c),  and  (e))  and  winter  wind  stress  ((b),  (d),  and  (f)). 
Positive  values  indicate  that  heat  is  transferred  from  the  ocean  to  the  atmosphere. 


To  explore  how  the  local  thermal  and  wind  forcing 
affect  the  interannual  variability  of  the  YSCWM,  the  inter¬ 
annual  components  of  winter  heat  flux  and  wind  stress  are 
decomposed  using  EOFs  (Figure  5).  EOFl  of  winter  heat  flux 
shows  a  positive  loading  over  the  total  Yellow  Sea,  which 
explains  68.6%  of  the  total  variance  (Figure  5(a)).  The  EOF2 
of  winter  heat  flux  represents  an  east- west  dipole  structure, 
with  positive  loading  in  the  east  and  negative  loading  in 
the  west  (Figure  5(c)).  The  EOF3  has  a  south-north  dipole 
structure  with  positive  loading  in  the  northern  Yellow  Sea 
and  negative  loading  in  the  southern  Yellow  Sea  (Figure  5(e)). 
EOFl  of  winter  wind  stress  is  characterized  by  negative 
values  over  almost  the  whole  Yellow  Sea,  which  account  for 
58.2%  of  the  total  variance  (Figure  5(b)).  The  EOF2  of  winter 
wind  stress  reveals  a  south-north  dipole  pattern  with  positive 
values  in  the  northern  Yellow  Sea  and  negative  values  in 
the  southern  Yellow  Sea  (Figure  5(d)).  The  EOF3  exhibits 
a  meridional  alternate  pattern  with  positive  loading  in  the 


Chinese  coast  and  negative  loading  in  the  Korean  coast  (Fig¬ 
ure  5(f)). 

The  EOFl  of  winter  heat  flux  is  negatively  (positive) 
correlated  with  the  A-YSCWM  (T-YSCWM)  (Figure  6(a)  and 
Table  2).  This  suggests  that  strong  winter  heat  loss  from  ocean 
to  atmosphere  in  the  Yellow  Sea  can  induce  large  A-YSCWM 
and  low  T-YSCWM.  EOFl  and  EOF3  of  winter  wind  stress 
evidently  affect  the  features  of  the  YSCWM  (Figure  6(c)  and 
Table  3).  It  is  shown  that  both  EOFl  and  EOF3  of  winter  wind 
stress  are  negatively  (positively)  correlated  with  A-YSCWM 
(T-YSCWM).  This  indicates  that  enhanced  winter  winds 
over  the  Yellow  Sea  tend  to  make  A-YSCWM  (T-YSCWM) 
increase  (decrease)  (Figure  5(a)). 

Figures  7(a)  and  7(c)  show  the  interannual  time  series 
of  winter  AO,  MEI,  and  PDO  indices,  as  well  as  the  EOFl, 
EOF2,  and  EOF3  of  heat  flux  and  wind  stress.  Only  the  winter 
AO  index  is  positively  correlated  with  the  EOF3  of  winter 
heat  flux  (Table  4).  The  variations  of  local  winter  heat  flux 
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Figure  6:  Interannual  time  series  of  (a)  EOFl,  EOF2,  and  EOF3  of  winter  heat  flux,  and  the  A-YSCWM  and  T-YSCWM  in  August  (b)  EOFl, 
EOF2,  and  EOF3  of  summer  heat  flux,  and  the  A-YSCWM  and  T-YSCWM  in  August  (c)  EOFl,  EOF2,  and  EOF3  of  winter  wind  stress,  and 
the  A-YSCWM  and  T-YSCWM  in  August  (d)  EOFl,  EOF2,  and  EOF3  of  summer  wind  stress,  and  the  A-YSCWM  and  T-YSCWM  in  August. 
All-time  series  are  standardized.  The  A-YSCWM  and  T-YSCWM  represent  the  area  and  averaged  temperature  of  the  YSCWM,  respectively. 


Table  1:  Correlation  coefficients  between  the  AO,  MEI,  PDO, 
and  the  YSCWM  area  and  temperature  in  August,  including  the 
interannual  and  decadal  components  in  both  summer  and  winter. 
Values  in  bold  are  statistically  significant  at  the  95%  confidence  level. 


Area 

Average  temperature 

Interannual  variability 

Winter  AO 

-0.23 

0.25 

Winter  MEI 

-0.27 

0.10 

Winter  PDO 

0.31 

-0.32 

Summer  AO 

-0.27 

0.30 

Summer  MEI 

-0.08 

0.14 

Summer  PDO 

-0.15 

0.12 

Decadal  variability 

Winter  AO 

-0.49 

-0.05 

Winter  MEI 

-0.41 

0.32 

Winter  PDO 

-0.43 

0.31 

Summer  AO 

-0.29 

-0.16 

Summer  MEI 

-0.61 

0.38 

Summer  PDO 

-0.26 

0.04 

and  wind  stress  are  less  related  to  the  winter  MEI  and  PDO 
(Table  4).  In  other  words,  the  winter  MEI  and  PDO  have 
less  influence  on  the  YSCWM  through  the  modulation  of  the 
winter  heat  flux  and  wind  stress. 

4.1.2.  Relationship  to  the  Climate  Forcings  in  Summer.  A 
positive  correlation  exists  between  the  summer  AO  index 
and  T-YSCWM,  with  a  correlation  coefficient  of  0.30  that 
is  statistically  significant  at  the  confidence  level  of  95%. 
The  correlations  between  MEI,  PDO  and  A-YSCWM,  T- 
YSCWM  are  low  and  not  statistically  significant  (Eigure  4(b) 
and  Table  1).  This  suggests  that  ENSO  and  PDO  may  not  be 
the  main  factors  controlling  the  interannual  variability  of  the 
YSCWM  in  summer. 

The  interannual  variability  of  summer  heat  flux  and  wind 
stress  is  analyzed  using  EOEs.  Eigures  6(b)  and  6(d)  show 
the  interannual  time  series  of  summer  AO,  MEI,  and  PDO 
indices,  as  well  as  the  EOEl,  EOP2,  and  EOP3  of  heat  flux  and 
wind  stress.  The  local  heat  flux  and  wind  stress  in  summer 
have  little  impact  on  the  variability  of  the  YSCWM  (Tables 
2  and  3).  Actually,  none  of  the  three  climate  indices  are 
correlated  with  the  local  heat  flux  and  wind  stress  (Eigures 
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Figure  7:  Interannual  time  series  of  (a)  winter  AO,  MEI,  PDO,  EOFl,  EOF2,  and  EOF3  of  winter  heat  flux;  (b)  summer  AO,  MEI,  PDO,  EOFl, 
EOF2,  and  EOF3  of  summer  heat  flux;  (c)  winter  AO,  ENSO,  PDO,  EOFl,  EOF2,  and  EOF3  of  winter  wind  stress;  (d)  summer  AO,  MEI,  and 
PDO,  and  EOFl,  EOF2,  and  EOF3  of  summer  wind  stress.  EOFl,  EOF2,  and  EOF3  of  heat  flux  and  wind  stress  are  standardized. 


7(b),  7(d)  and  Table  5).  This  means  that  summer  AO,  ENSO, 
and  PDO  do  not  affect  the  YSCWM  through  modulating  the 
local  heat  flux  and  wind  stress. 

4.2.  Decadal  Variability 

4.2.1.  Relationship  with  the  Climate  Forcings  in  Winter.  The 
decadal  variability  of  the  A- YSCWM  has  a  high  correlation 
with  the  winter  AO  index.  The  negative  correlation  between 
the  winter  AO  and  the  area  of  the  YSCWM  on  the  decadal 
time  scale  is  about  -0.49  (Figure  8(a)  and  Table  1).  When 
the  winter  AO  index  is  in  its  negative  phase  during  1958- 
1988,  A-YSCWM  usually  has  positive  anomaly.  After  1988, 
the  AO  index  is  continuously  positive,  and  A-YSCWM  shows 
a  persistent  negative  anomaly.  When  the  AO  index  is  in  its 
positive  phase,  the  heat  flux  loading  associated  with  EOFl  is 
negative,  indicating  a  net  heat  flux  from  the  atmosphere  and 
into  the  ocean.  A-YSCWM  also  reduces  accordingly  (Figures 
9, 11(a)  and  Table  4). 

The  PDO  index  is  correlated  with  A-YSCWM  and  T- 
YSCWM,  with  the  correlation  coefficient  of  -0.43  and  0.31, 
respectively  (Figure  8(a)  and  Table  1),  and  maintains  positive 


or  negative  values  prevailing  for  20-30  year  periods  [21,  22]. 
The  decadal  cycle  of  the  YSCWM  also  shows  a  similar  feature 
(Figure  8(a)).  In  winter,  the  PDO  index  has  strong  links 
with  the  Aleutian  low.  Negative  phase  of  the  PDO  represents 
warm  anomalies  in  the  extratropical  North  Pacific  [20].  The 
spatial  pattern  of  the  PDO  is  strongly  associated  with  the 
EOFl  of  heat  flux  (Figure  10(a)).  When  the  PDO  is  in  its 
positive  phase,  the  heat  flux  from  atmosphere  to  ocean  is 
increased  in  the  Yellow  Sea,  which  reduces  the  A-YSCWM 
(Figures  9, 11(a)  and  Table  4).  The  reverse  is  true  during  the 
negative  PDO  phase.  The  PDO  also  has  a  connection  with 
EOF2  and  EOF3  of  the  wind  stress  (Figure  10(c)).  Positive 
phase  of  the  PDO  induces  positive  wind  stress  anomaly  in  the 
Yellow  Sea  (Figures  9, 10(c)  and  Table  4).  Due  to  the  shallow 
depth  of  the  Yellow  Sea,  a  stronger  wind  stress  can  strengthen 
mixing,  which  produces  a  smaller  A-YSCWM  (Figure  11(c) 
and  Table  3). 

The  correlations  between  MEI  and  A-YSCWM  and  T- 
YSCWM  are  -0.41  and  0.32,  respectively  (Eigure  8(a)  and 
Table  1).  From  1976  to  2006,  the  MEI  represents  the  warm 
ENSO  phase.  The  corresponding  T-YSCWM  (A-YSCWM) 
is  a  positive  (negative)  anomaly.  However,  the  correlations 
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Figure  8:  Decadal  time  series  of  (a)  winter  AO,  MEI,  and  PDO 
index,  and  the  A-YSCWM  and  T-YSCWM  in  August  (b)  summer 
AO,  MEI,  and  PDO  index,  and  the  A-YSCWM  and  T-YSCWM 
in  August.  Time  series  of  the  A-YSCWM  and  T-YSCWM  are 
standardized.  The  A-YSCWM  and  T-YSCWM  represent  the  area 
and  averaged  temperature  of  the  YSCWM,  respectively. 

Table  2:  Correlation  coefficients  between  EOFl,  EOF2,  and  EOF3 
of  heat  flux  and  A-YSCWM,  T-YSCWM  in  August,  including  the 
interannual  and  decadal  components  in  both  summer  and  winter. 
Bold  texts  indicate  that  the  correlations  are  statistically  significant  at 
the  95%  confidence  level. 


A-YSCWM 

T-YSCWM 

Interannual  variability 

EOFl  of  winter  heat  flux 

-0.32 

0.37 

EOF2  of  winter  heat  flux 

-0.16 

0.14 

EOF3  of  winter  heat  flux 

0.04 

-0.08 

EOFl  of  summer  heat  flux 

-0.13 

0.07 

EOF2  of  summer  heat  flux 

0.04 

0.03 

EOF3  of  summer  heat  flux 

0.17 

0.01 

Decadal  variability 

EOFl  of  winter  heat  flux 

0.74 

-0.08 

EOF2  of  winter  heat  flux 

0.02 

-0.33 

EOF3  of  winter  heat  flux 

0.02 

-0.39 

EOFl  of  summer  heat  flux 

0.59 

-0.64 

EOF2  of  summer  heat  flux 

-0.56 

-0.11 

EOF3  of  summer  heat  flux 

0.22 

-0.12 

between  MEI  and  local  wind  stress  and  heat  flux  are  low 
(Figures  10(a),  10(c)  and  Table  4).  The  results  indicate  that 
the  winter  ENSO  event  is  strongly  linked  to  the  decadal 
variability  of  the  YSCWM,  but  not  through  changes  in  wind 
stress  and  heat  fluxes. 


Table  3:  Correlation  coefficients  between  EOFl,  EOF2,  and  EOF3 
of  wind  stress  and  A-YSCWM,  T-YSCWM  in  August,  including  the 
interannual  and  decadal  components  in  both  summer  and  winter. 
Bold  texts  indicate  that  the  correlations  are  statistically  significant  at 
the  95%  confidence  level. 


A-YSCWM 

T-YSCWM 

Interannual  variability 

EOFl  of  winter  wind  stress  -0.33 

0.33 

EOF2  of  winter  wind  stress 

0.08 

-0.09 

EOF3  of  winter  wind  stress 

-0.31 

0.27 

EOFl  of  summer  wind  stress 

-0.07 

0.11 

EOF2  of  summer  wind  stress 

0.15 

-0.07 

EOF3  of  summer  wind  stress 

-0.08 

0.12 

Decadal  variability 

EOFl  of  winter  wind  stress  -0.13 

-0.23 

EOF2  of  winter  wind  stress 

-0.56 

0.75 

EOF3  of  winter  wind  stress 

0.31 

0.08 

EOFl  of  summer  wind  stress 

0.19 

0.56 

EOF2  of  summer  wind  stress 

-0.19 

0.14 

EOF3  of  summer  wind  stress 

-0.39 

0.48 

4.2.2.  Relationship  to  Climate  Forcings  in  Summer.  On  the 
decadal  time  scale,  the  A-YSCWM  is  negatively  correlated 
with  the  summer  AO  index,  with  a  correlation  coefficient 
of  -0.30  (Figure  8(b)  and  Table  1).  The  A-YSCWM  and  the 
summer  MEI  are  negatively  correlated  with  a  coefficient 
of  -0.60  (Figure  8(b)  and  Table  1).  The  correlation  between 
T-YSCWM  and  summer  MEI  is  positive  with  the  corre¬ 
lation  coefficient  of  0.38  (Figure  8(b)  and  Table  1).  On  a 
decadal  time  scale,  the  summer  MEI  during  1966-1978  is 
in  its  negative  phase  while  the  anomaly  in  A-YSCWM  (T- 
YSCWM)  is  positive  (negative).  After  1978,  the  summer 
MEI  is  in  its  positive  phase  (Figure  8(b)).  Correspondingly, 
the  A-YSCWM  (T-YSCWM)  anomaly  is  negative  (positive). 
However,  when  PDO  turns  to  a  weak  signal  in  summer,  the 
correlation  between  the  PDO  and  the  YSCWM  becomes  not 
significant  (Figure  8(b)  and  Table  1). 

The  summer  AO  and  MEI  may  affect  both  EOF2  and 
EOF3  of  the  summer  heat  flux  and  EOFl  of  summer  wind 
stress,  modulating  the  decadal  variability  of  the  YSCWM 
(Figures  10(b),  10(d)  and  Table  5).  However,  because  of  the 
low  correlation  between  the  EOF2  of  the  summer  heat  flux 
and  the  decadal  variability  of  the  YSCWM,  AO  and  MEI 
may  not  affect  the  decadal  variability  of  the  YSCWM  through 
changes  in  the  EOF2  of  the  summer  heat  flux  (Table  3).  The 
summer  AO  and  ENSO  indices  are  negatively  correlated  with 
EOF3  of  summer  heat  flux  and  EOFl  of  summer  wind  stress 
(Tables).  When  the  summer  AO  and  ENSO  are  in  their 
positive  phases,  the  heat  flux  from  the  atmosphere  into  the 
ocean  is  enhanced  (Figures  9,  10  and  Table  5).  Therefore, 
the  A-YSCWM  decreases,  and  the  T-YSCWM  increases 
(Figure  11  and  Table  3). 

5.  Summary 

We  use  the  CORA  reanalysis  data  to  investigate  the  low- 
frequency  variability  of  the  YSCWM.  First,  we  examine  the 
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Figure  9:  EOFl,  EOF2,  and  EOF3  of  decadal  components  for  winter  heat  flux  ((a),  (e),  and  (i))  and  summer  heat  flux  ((b),  (f),  and  (g)),  winter 
wind  stress  ((c),  (g),  and  (k)),  and  summer  wind  stress  ((d),  (h),  and  (1)).  Positive  values  indicate  that  heat  is  transferred  from  the  ocean  to 
the  atmosphere. 


Table  4:  Correlation  coefficients  between  EOFl,  EOF2,  and  EOF3  of  winter  heat  flux,  wind  stress,  and  winter  AO,  ENSO,  and  PDO,  including 
the  interannual  and  decadal  components.  Bold  texts  indicate  that  the  correlations  are  statistically  significant  at  the  95%  confidence  level. 


Winter  AO 

Winter  MEI 

Winter  PDO 

EOFl  of  winter  heat  flux 

0.26 

0.02 

-0.11 

EOF2  of  winter  heat  flux 

-0.08 

0.10 

-0.11 

EOF3  of  winter  heat  flux 

0.47 

0 

-0.03 

Interannuai  variability 

EOFl  of  winter  wind  stress 

0.09 

0.09 

-0.01 

EOF2  of  winter  wind  stress 

0.10 

-0.05 

-0.04 

EOF3  of  winter  wind  stress 

0.04 

0.17 

-0.11 

EOFl  of  winter  heat  flux 

-0.59 

-0.12 

-0.47 

EOF2  of  winter  heat  flux 

-0.09 

-0.06 

0.13 

EOF3  of  winter  heat  flux 

0.57 

0.08 

-0.31 

Decadal  variability 

EOFl  of  winter  wind  stress 

0.17 

-0.04 

-0.10 

EOF2  of  winter  wind  stress 

-0.21 

0.25 

0.33 

EOF3  of  winter  wind  stress 

-0.14 

-0.12 

-0.77 

formation  and  evolution  of  the  YSCWM  on  the  seasonal 
scale.  The  YSCWM  forms  in  spring  and  reaches  its  maximum 
extension  in  summer.  Since  the  early  fall,  the  area  occupied  by 
the  YSCWM  starts  to  decrease.  In  December,  the  YSCWM 
disappears.  The  evolution  process  of  the  YSCWM  derived 
from  CORA  is  consistent  with  the  evolution  showed  by  the 
observations  as  described  in  previous  studies  [24]. 


Next,  we  examine  the  interannual  and  decadal  variability 
of  the  YSCWM  and  its  association  with  the  climate  forcing 
including  AO,  ENSO,  and  PDO.  The  climate  forcing  plays  an 
important  role  in  modulating  the  YSCWM  variability.  For 
the  interannual  variability  of  the  YSCWM,  the  correlation 
between  the  winter  PDO  and  the  YSCWM  is  the  strongest 
among  the  three  climate  indices,  which  indicates  that  the 
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Figure  10:  Decadal  time  series  of  (a)  winter  AO,  MEI,  PDO,  EOFl,  EOF2,  and  EOF3  of  winter  heat  flux;  (b)  summer  AO,  MEI,  PDO,  EOFl, 
EOF2,  and  EOF3  of  heat  flux;  (c)  winter  AO,  MEI,  PDO,  and  EOFl,  EOF2,  and  EOF3  of  winter  wind  stress;  (d)  summer  AO,  MEI,  PDO,  and 
EOFl,  and  EOF2  and  EOF3  of  summer  wind  stress.  EOFl,  EOF2,  and  EOF3  of  heat  flux  and  wind  stress  are  standardized. 


Table  5:  Correlation  coefficients  between  EOFl,  EOF2,  and  EOF3  of  summer  heat  flux,  wind  stress,  and  winter  AO,  ENSO,  and  PDO,  including 
the  interannual  and  decadal  components.  Bold  texts  indicate  that  the  correlations  are  statistically  significant  at  the  95%  confidence  level. 


Summer  AO 


Summer  MEI 


Summer  PDO 


Interannual  variability 


EOFl  of  summer  heat  flux 
EOF2  of  summer  heat  flux 
EOF3  of  summer  heat  flux 
EOFl  of  summer  wind  stress 
EOF2  of  summer  wind  stress 
EOF3  of  summer  wind  stress 


-0.07 

0.04 

-0.01 

0.01 

-0.04 

-0.01 


0.01 

0.09 

0.17 

0.01 

0.11 

0.11 


0.01 

0.13 

0.12 

-0.05 

0.12 

0.15 


Decadal  variability 


EOFl  of  summer  heat  flux 
EOF2  of  summer  heat  flux 
EOF3  of  summer  heat  flux 
EOFl  of  summer  wind  stress 
EOF2  of  summer  wind  stress 
EOF3  of  summer  wind  stress 


0.08 

0.41 

-0.40 

-0.57 

0.15 

-0.02 


-0.13 

0.47 

-0.61 

-0.36 

-0.07 

0.31 


0.14 

0.26 

-0.79 

-0.48 

-0.03 

0.04 


effect  of  the  winter  PDO  on  the  YSCWM  is  dominant.  In 
this  study,  we  do  not  investigate  how  PDO  modulate  the 
Yellow  Sea  circulation  to  influence  the  YSCWM.  Previous 
studies  suggested  that  the  signal  of  PDO  can  be  transferred 
into  the  Yellow  Sea  through  Kuroshio  and  therefore  affect  the 


circulation  of  the  Yellow  Sea  [16,  25].  It  was  also  suggested 
that  Yellow  Sea  warm  current  can  influence  the  variability  of 
the  YSCMW  [26].  It  is  necessary  to  further  investigate  in  a 
future  study  how  PDO  influences  the  Yellow  Sea  circulation 
and  modulates  the  YSCWM.  It  is  also  found  that  the  local 
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Figure  11:  Decadal  time  series  of  (a)  EOFl,  EOF2,  and  EOF3  of  winter  heat  flux,  and  the  A-YSCWM  and  T-YSCWM  in  August  (b)  EOFl, 
EOF2,  and  EOF3  of  summer  heat  flux,  and  the  A-YSCWM  and  T-YSCWM  in  August  (c)  EOFl,  EOF2,  and  EOF3  of  winter  wind  stress,  and  the 
A-YSCWM  and  T-YSCWM  in  August  (d)  EOFl,  EOF2,  and  EOF3  of  summer  wind  stress,  and  the  A-YSCWM  and  T-YSCWM  in  August.  All 
of  the  time  series  are  standardized.  The  A-YSCWM  and  T-YSCWM  represent  the  area  and  averaged  temperature  of  the  YSCWM,  respectively. 


wind  stress  and  heat  flux  in  summer  have  little  influence  on 
the  variability  of  the  YSCWM.  In  summer,  no  correlation  is 
found  between  the  climate  indices  and  the  local  heat  flux  and 
wind  stress.  This  suggests  that  the  AO,  ENSO,  and  PDO  may 
not  affect  the  inter  annual  variability  of  the  YSCWM  through 
modulating  the  local  heat  flux  and  wind  stress  in  summer. 

For  the  decadal  variability  of  the  YSCWM,  the  negative 
correlations  between  the  winter  AO,  PDO,  and  the  area  of  the 
YSCWM  are  significant.  The  winter  AO  mainly  controls  the 
first  mode  of  heat  flux  to  modulate  the  decadal  variability  of 
the  YSCWM.  When  the  winter  AO  is  in  its  positive  phase, 
the  heat  flux  from  atmosphere  to  ocean  in  the  northwestern 
Yellow  Sea  is  abnormally  high.  The  average  temperature  of 
the  YSCWM  also  increases  accordingly.  The  winter  PDO  is 
strongly  connected  to  the  first  mode  of  winter  heat  flux.  When 
the  winter  PDO  is  in  its  positive  phase,  the  increase  of  heat 
flux  from  atmosphere  to  ocean  in  the  Yellow  Sea  can  reduce 
the  area  of  the  YSCWM.  At  the  same  time,  the  winter  PDO 
has  positive  correlation  with  the  sea  surface  pressure.  Along 
with  the  PDO  enhanced,  the  enhancement  of  sea  surface 
pressure  is  in  conformity  with  wind  speed.  The  positive  PDO 
can  induce  positive  wind  stress  anomaly  in  the  Yellow  Sea 
through  the  modulation  of  the  EOF2  and  EOF3.  Due  to  the 


shallow  depth  of  the  Yellow  Sea,  the  strong  wind  stress  and  the 
enhanced  sea  surface  pressure  can  strengthen  mixing,  reduc¬ 
ing  the  area  of  the  YSCWM.  The  correlation  between  winter 
MEI  and  the  YSCWM  is  significant,  suggesting  a  link  with 
decadal  variability  of  the  YSCWM.  However,  this  relationship 
cannot  be  explained  by  changes  in  the  wind  stress  and  heat 
fluxes.  In  summer,  ENSO  has  the  strongest  influence  on  the 
YSCWM  variability.  ENSO  may  exert  influence  on  the  third 
mode  of  summer  heat  flux  and  the  first  mode  of  summer  wind 
stress  to  modulate  the  decadal  variability  of  the  YSCWM. 
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